System Size, Energy and Centrality Dependence of Pseudo-rapidity Distributions of 
Charged Particles in Relativistic Heavy Ion Collisions 
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We present the first measurements of the pseudorapidity distribution of primary charged particles 
in Cu+Cu collisions as a function of collision centrality and energy, ys NN = 22.4, 62.4 and 200 GeV, 
over a wide range of pseudorapidity, using the PHOBOS detector. Making a global comparison of 
Cu + Cu and Au + Au results, we find that the total number of produced charged particles and the 
rough shape (height and width) of the pseudorapidity distributions are determined by the number 
of nucleon participants. More detailed studies reveal that a more precise matching of the shape of 
the Cu + Cu and Au + Au pseudorapidity distributions over the full range of pseudorapidity occurs 
for the same N par t/2A value rather than the same N pa tt value. In other words, it is the collision 
geometry rather than just the number of nucleon participants that drives the detailed shape of the 
pseudorapidity distribution and its centrality dependence at RHIC energies. 

PACS numbers: 25.75.-q, 25.75.Dw 



The advent of Cu + Cu collisions from the Relativistic 
Heavy Ion Collider (RHIC) at energies similar to those of 
the earlier Au + Au collisions presents a new opportunity 
to measure the system size dependence of important 
observables using different collision geometries. The 
Cu + Cu results are expected to provide critical tests of 
the parametric dependence of the pseudorapidity density 
of charged particles, dN c h/d7y, observed previously in 
Au + Au collisions p], 0, HJ . They significantly extend 
the range of measurements with respect to the number 
of participant nucleons, N part , compared to Au + Au and 
also allow for a direct comparison at the same N part . 

The observed dN c h/d?7 is a conceptually well-defined 
quantity that reflects most effects that contribute to 
particle production in heavy- ion collisions. It is sensitive 
to the initial conditions of the system, i.e. parton shad- 
owing, and also to the effects of rescattering and hadronic 
final-state interactions. In short, the full distribution of 
dN c h/dry represents a time- integral of particle production 
throughout the entire heavy-ion collision. 

In this letter, we present the first measurements of the 
dN c h/d?7 of primary charged particles over a broad range, 
1 77 1 < 5.4, for Cu + Cu collisions at a variety of collision 



centralities. These measurements were made in the same 
detector for s /s^= 22.4, 62.4 and 200 GeV allowing for 
a reliable systematic study of particle production as a 
function of energy. The Cu + Cu results are compared 
to results from Au + Au 0, 0] and d+Au collisions [1, 
5] at similar energies obtained with the same detector. 
This led us to perform a comprehensive examination of 
particle production in Cu + Cu and Au + Au collisions 
for the same number of nucleon participant pairs and for 
the same fraction of total cross section in both systems, 
and to study the interplay between energy and system 
size. 

The data were obtained with the multiplicity array 
of the PHOBOS detector @] during the RHIC 2005 
run. This array is identical to that used in our study 
of Au + Au and d+ Au collisions 0, S 0, H . Monte 
Carlo simulations of the detector performance were based 
on the HIJING event generator @] and GEANT [| 
simulations, folding in the signal response for scintillator 
counters (Paddles, covering 3.0 < \r]\ < 4.5), and silicon 
sensors. 



The Au + Au data at 



19.6, 62.4 and 200 GeV 



used in this paper for comparison with Cu + Cu results 
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FIG. 1: dN c h/d»7 distributions of primary charged particles 
from Cu + Cu collisions at ^/s^7 = 22.4, 62.4, and 200 GeV 
for the specified centrality bins. The typical systematic errors 
(90% C.L.) are shown as bands for selected centrality bins. 



FIG. 2: dN ch /d?7 distributions in Cu+Cu and Au+Au @, 0] 
collisions selected to yield similar (N par t), a) at 200 GeV, b) 
at 62.4 GeV and c) at 22.4 GeV. The grey band indicates 
the systematic uncertainty (90% C.L.) for Cu+Cu. Errors for 
Au + Au are not shown for clarity. 



were taken with the PHOBOS detector @, E|. In the 
present work, the Au + Au data at 19.6 GeV have been 
reanalyzed using an improved treatment of dead detector 
channels, resulting in slightly smaller multiplicity than 
given in previous publications. However, the dN c h/d?y 
distributions shown in the present work and those pub- 
lished in Refs. 0, Q agree within the systematic errors. 
The 19.6 Au + Au data and the new Cu + Cu data were 
analyzed in the same way as the previously published 
data of Au + Au at 200 and 62.4 GeV @, If, using two 
analysis methods @, "hit-counting" and "analog". In 
the present work, the analog results were corrected for 
the over-counting of multiply-charged fragments such as 
He emitted in the forward regions. The correction, which 
was largest for the lowest energies and most peripheral 
collisions, changed the total number of charged particles 
by less than 6%, and has been taken into account in the 
systematic error assignments. 

The centrality determination procedure applied for 
Cu + Cu collisions is the same as for Au + Au collisions 
at similar energies @, i]. For the 200 and 62.4 GeV 
data sets, the centrality was estimated from the data 



using the truncated mean of the Paddles signals. Using 
several methods, based on the HIJING [7] and AMPT 
[Io| models, we have estimated our minimum-bias trigger 
efficiency for events with a vertex near the nominal inter- 
action point to be 84 ± 5% and 75 ± 5% in Cu + Cu col- 
lisions at 200 and 62.4 GeV, respectively. At the lowest 
energies, 22.4 GeV (Cu + Cu) and 19.6 GeV (Au + Au), 
the much lower beam rapidity (ybeam ~ 3) precludes the 
use of the Paddles for centrality determination. Instead, 
we construct a different quantity, "EOCT" , the path- 
length-corrected sum of the energy deposited in the Oc- 
tagon (silicon) detector (|ry| < 3.2V This procedure has 
been discussed in detail in Ref. 0. We use a Glauber 
model calculation implemented in HIJING simulations 
to estimate (N part ) for each centrality bin. The (N part ) 
values for various centrality bins for Cu + Cu collisions 
are given in Table |TJ 

Figure [T] shows the primary charged particle dN c h/d?7 
distributions measured in Cu + Cu collisions at ^/s NN = 
22.4, 62.4, and 200 GeV for different centrality bins. 
The statistical errors are negligible. Both the height and 
width of the dN c h/d?7 distributions increase as a function 
of energy as has been seen for Au + Au collisions • 
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3: Cu + Cu and Au + Au U data at RHIC energies, 




Cu+Cu Au+Au 

0-3% 0-10% 

200 GeV 0.84 0.83 

62.4 GeV 0.81 0.80 
22.4 GeV 0.81 



plotted as dN ch /drj' /{N part /2), where rf = r\ - y beam , for 
a) 0-6% most central events and b) events with similar 
value of N par t/2A. For the 0-6% central collisions of Cu + Cu 
(Au + Au) at 200, 62.4 and 22.4 (19.6) GeV, the values of 
N P art/2A are 0.82 (0.87),0.78 (0.85) and 0.78 (0.86), respec- 
tively. Systematic errors (90% C.L.) are shown for typical 
points. 



The comparison of dN c h/dr? distributions for Cu + Cu 
and Au + Au collisions at the same energy, for centrality 
bins chosen so that the average number of participants 
in both systems is similar, is presented in Fig. [2j No 
scaling factors are applied. At 200 GeV, Fig. [2^, the 
dN c h/dr/ distributions in both systems at similar (N part ) 
agree within systematic errors, both in height and width. 
At 62.4 and 22.4 GeV, Fig. [2p and c, we observe that 
the distributions agree within systematic errors at midra- 
pidity but not in the fragmentation regions (i.e. high 
|ry|). The dN crl /d?7 distributions of Au + Au collisions 
at 19.6 and 62.4 GeV have been interpolated linearly 
in ln(^/s NN ) to obtain the scaled Au + Au data at 22.4 
GeV. The fragmentation region dependence on the size 
of the colliding nuclei is exhibited through increased 
charged particle production in Au + Au as compared to 
Cu + Cu collisions which may be attributed to the two 
excited nuclear remnants being bigger in Au + Au than 
in Cu + Cu collisions. This effect is most visible at the 
lowest energies where the broad r\ coverage gives access 
to \rj\ > Vb earn, - 

The interplay between the system size and the collision 
energy is shown in Fig. [3] It shows the scaled pseu- 
dorapidity particle densities, dN crl /d?//(N par t/2) (where 
77' = 77 — ybeam corresponds effectively to the rest frame 
of one of the colliding nuclei [HQ) for Cu + Cu and 
Au + Au collisions for centrality bins with a) the same 
fraction of total cross section 0-6% and b) similar value 
of N part /2A, where A is the mass number of the colliding 
nuclei. It should be noted that systems with matching 

part 



Npart /2A values will also have matching N spec /N 
values, where N spec = 2 A - N part is the number of non- 
participating nucleons (spectators). The results shown 
in Fig. [3] suggest that in symmetric nucleus- nucleus col- 
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FIG. 4: The ratio, Rpg art (V), of dN ch /drf / '(N part /2) for 
Cu + Cu and Au + Au collisions at RHIC energies, a) com- 
paring the 35-40% bin to the 0-6% most central bin and 
b) Rpc art ('/) f° r Cu + Cu and Au + Au for centrality bins 
selected such that N part /2A is similar for the two systems. 

The inset figure (a) represents the RpQ art (»/) only for Au + Au 
data for four different energies fj|. The errors represent a 90% 
C.L. systematic error on the ratio. Note : the scaled Au + Au 
at 22.4 GeV has not been added in the figure because the 
ratio, Rpg art (7/), for a given system is independent of energy. 
The values of N par t/2A are given in Fig. [3] and the value 
of Np ar t/2 A for Cu + Cu and Au + Au for 35-40% collisions 
centrality at 22.4 (19.6) GeV are both equal to 0.24. 



lisions the particle density per nucleon participant pair 
at the midrapidity region does not depend on the size 
of the two colliding nuclei but only on the collision 
energy and geometry. In the fragmentation region, the 
phenomenon of extended longitudinal scaling observed 
in Au + Au 0, Q and d + Au [H collisions is also 
present in the Cu + Cu data. The Cu + Cu and Au + Au 
collisions exhibit the same extended longitudinal scaling 
in central collisions. This suggests that the extended 
longitudinal scaling holds independent of the collision 
energy and system size. The comparison presented in 
Fig. [3] reveals an interesting feature: for centrality bins 
corresponding to the same fraction of total cross section 
in Cu + Cu and Au + Au collisions and at similar energy, 
the dN c h/d?7/ (N part /2) distributions in both systems 
agree within errors over the full range of rj. Slightly 
better agreement of dN cri /d?7/(N par t/2) distributions is 
obtained for centrality bins selected to yield similar value 
of N part /2A in both systems as presented in Fig. [3)3. 

To study the centrality dependence of extended lon- 
gitudinal scaling in Cu + Cu and Au + Au collisions, we 
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FIG. 5: The integrated total of primary charged particle 
production obtained by extrapolating the data at each energy 
into the unmeasured region is shown as a function of centrality 
in Cu + Cu collisions. The Au + Au, d+Au and inelastic 
p(p) + p data are taken from Refs. 3, 5]. The uncertainty on 

part 



N^' and N par t has been included in the error bars 



examine the ratios of dN c h/aV/(Np art /2) for central and 

semi-central, denoted Rp£T rt , at different energies as a 
function of rf . The inset of Fig. [4ji shows the previously 
published results for Au + Au collisions [3j , indicating 
that the change in shape as a function of centrality is in- 
dependent of beam energy. The ratios Rpc* rt for Cu + Cu 
at three energies exhibit the same feature (solid points) 
in Fig. [4Jl. By comparing to Au + Au results (open 
points) this ratio is found to be similar in the midrapidity 
region —2.5 <rf< 0.5 but in the region r/ > 0.5 the 
ratio for Au + Au is higher. To study this difference, 
we plot in Fig. HJd the Rpg"' ratio for Cu + Cu and 
Au + Au for centrality bins selected to represent similar 
initial geometry even more precisely, i.e. similar value 
of N part /2A for the two systems. Using this comparison 
criterion we observe good agreement between the two 
systems over the full range of pseudorapidity. 

The values of total charged particle multiplicity, Nj^ * , 
estimated by extrapolation to the unmeasured region of 
pseudorapidity, are presented for Cu + Cu collisions at 
v /s^7 = 22.4, 62.4 and 200 GeV in Fig. [5] and given in 
Tabled] The method used to determine N^* is detailed 
in Ref. [ll| for Au + Au collisions. At all energies, N^' 
is obtained by averaging the results of two techniques. 
One involved fitting a Wood-Saxon functional form to 
the data for 1 77 1 < 8 at the two higher energies and 
M < yboam at the lowest energy where y£ e u ^ (22.4 GeV) 
= 3.2 and y£^{19.Q GeV) = 3.0. The other involved 
simply integrating the lowest energy data and using the 



extended longitudinal scaling result to extrapolate the 
higher energy data into the unmeasured regions. The 
Cu + Cu and Au + Au results are compared at the same 
energies, 62.4 and 200 GeV, as well as at nearly the same 
energy 22.4 (CuCu) and 19.6 GeV (Au + Au) in Fig. 
We observe that N^' scales approximately linearly with 
(Np ar t) in both Cu + Cu and Au + Au collisions, and 
has similar values for the same (N part ). The comparison 
indicates that the transition between inelastic p(p) + p 
and Cu + Cu collisions is not controlled simply by the 
number of participants, as even the very central d + Au 
multiplicity per participant pair shows little sign of con- 
tinuity to the Cu + Cu results. 



In summary, the measured pseudorapidity 
distributions of charged particles and the estimated 
total charged particle multiplicity in Cu + Cu collisions 
are presented as a function of collision centrality and 
energy, ^/s^= 22.4, 62.4 and 200 GeV. The results show 
that Np ar t is the scaling variable unifying the centrality 
dependence for Cu + Cu and Au + Au in the midrapidity 
region and for N^*. However, the best agreement of the 
pseudorapidity distributions per nucleon participants 
over the full range of r\ between Cu + Cu and Au + Au 
in the central collisions at the same energy is obtained 
for centrality bins selected to yield similar value of 
Npart /2A in both systems. The Cu + Cu and Au + Au 
results at similar energy show that the particle density 
per nucleon participant pair in the midrapidity region is 
similar in both systems. The phenomenon of extended 
longitudinal scaling in Cu + Cu and Au + Au collisions 
holds independent of colliding energy and system size. A 
dependence on the size of the colliding nuclei is observed 
in the pseudorapidity distributions in the fragmentation 
region at low energies, 22.4 (Au + Au 19.6) and 62.4 
GeV, when the collision centrality of the two systems is 
selected for similar N par t • This may be attributed to the 
two excited nuclear remnants being bigger in Au + Au 
than in Cu + Cu collisions. The essential role of collision 
geometry when comparing pseudorapidity distributions 
of charged particles between nuclear species is clearly 
demonstrated. 
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TABLE I: The estimated number of nucleon participants, (N par t), and the total charged particle multiplicity, N^*, extrapolated 
to the unmeasured region for Cu+Cu collisions in different centrality bins are presented. All errors are systematic (90% C.L.). 



Centrality 


200 GeV 


62.4 GeV 


22.4 GeV 


Bin (%) 


(Np ar t) 


AT Tot 

N c h 


(N 


3 art) 


i\rTot 

N c h 


(N 


3 art) 


N^ 


fl— ^ 


106 ± 3 


1541 ± 70 


102 ± 3 


833 ± 36 


103 ± 3 


535 ± 


23 


o u 


100 ± 3 


1407 ± 68 


95 


± 


3 


781 ± 34 


95 


± 


3 


482 ± 


21 


U J-U 


91 ± 3 


1262 ± 59 


88 


± 


3 


721 ± 32 


86 


± 


3 


431 ± 


19 


10-15 


79 ± 3 


1084 ± 51 


76 


± 


3 


635 ± 27 


74 


± 


3 


375 ± 


18 


15-20 


67 ± 3 


917 ± 43 


65 


± 


3 


541 ± 24 


63 


± 


3 


320 ± 


15 


20-25 


57 ± 3 


771 ± 38 


55 


± 


3 


460 ± 21 


53 


± 


3 


273 ± 


14 


25-30 


47 ± 3 


645 ± 32 


47 


± 


3 


386 ± 17 


44 


± 


3 


230 ± 


13 


30-35 


40 ± 3 


538 ± 27 


38 


± 


3 


323 ±15 


37 


± 


3 


194 ± 


12 


35-40 


33 ± 3 


444 ± 23 


32 


± 


3 


270 ± 13 


30 


± 


3 


162 ± 


12 


40-45 


27 ± 3 


364 ± 19 


25 


± 


3 


223 ± 11 


25 


± 


3 


135 ± 


11 


45-50 


22 ± 3 


293 ± 15 


21 


± 


3 


183 ± 9 


20 


± 


3 


112 ± 


11 


50-55 


17 ± 3 


234 ± 13 


16 


± 


3 


147 ± 8 


16 


± 


3 


92 ± 


11 



